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Abstract 

The 2D flow of a foam confined in a Hele-Shaw cell through a con- 
traction is investigated. Its rheological features are quantified using 
image analysis, with measurements of the elastic stress, rate of plas- 
ticity, and velocity. The behavior of the velocity strongly differs at 
the contraction entrance, where the flow is purely convergent, and at 
the contraction exit, where a velocity undershoot and a re- focussing 
of the streamlines are unraveled. The yielded region, characterized by 
a significant rate of plasticity and a maximal stress amplitude, is con- 
centrated close to the contraction. These qualitative generic trends do 
not vary significantly with the flow rate, bubble area and contraction 
geometry, which is characteristic of a robust quasistatic regime. Using 
surfactants with a high surface viscoelasticity, a marked dependence 
of the elastic stress on the velocity is exhibited. The results show that 
the rate of plasticity does not only depend on the local magnitude of 
the deformation rate, but also crucially on the orientation of both elas- 
tic stresses and deformation rate. It is also discussed how the viscous 
friction controls the departure from the quasistatic regime. 



*Institut de Physique de Rennes, UMR 6251 CNRS/Universite de Rennes 
1, Batiment 11 A, Campus Beaulieu, 35042 Rennes Cedex, France. E-mail: 
benjamin.dollet@univ-rennesl .fr 



1 



1 Introduction 



Liquid foams are used in many industrial and domestic applications, such as 
ore notation, enhanced oil recovery or personal car e. Many of these uses take 



advan tage of the rich mechanical behavior of foams fjHohler and Cohen-Addad 



20051 ): under low applied strain or stress, they are elastic solids, whereas un- 



der high strain or stress, they undergo plastic flow, resulting of m any elemen- 
tary plastic events, the so-called Tls fjWeaire and Hutzlerl . Il999l ). character- 
ized by the topological rearrangement of four neighboring bubbles. Therefore, 
foams belong to the wide class of the complex fluids. Since their constitu- 
tive items, the bubbles, are of convenient size (typically 10~ 4 to 10 -2 m) 
for observation, they are particularly suited to relate a macroscopic mechan- 
ical response, measured for instance by rheometry, to the microstructural 
behavior, which is of paramount importance to inspire or test constitutive 
rheological models. 



Since bubbles are very efficient light scatterers, exper iments on 3D 



require techniques such as Diffusive- Wave Spectroscopy ( 


Durian et al. 


1991; 


GoDal and 


Durian. 1995: Earnshaw and Wilson 




1995 




Holder et al.. 


1997: 


Vera et al.. 


2001 


) or X-ray tomography ( 


Lambert et al.. 


2007 


). The former 



oams 



technique does not give precisely the bubble shape, and the latter is up 
to now limited by its long acquisition time. To overcome these difficulties, 
many experiment s have been p e rform ed on bubble monolayers, the so-called 
quasi-2D foams fjVaz and Coxl . 120051 ). where all bubbles c an be easily im- 



aged up to high frame rates. Since the seminal study of iDebregeas et al. 



(1200 lh . many fea tures of the flow of qua si-2D foams have b een quantified 



including elastic (lAsipauskas et all 120031 ) and pl astic aspects (Dennin 



Dollet and Graner. 



2007^ . and viscous stresses (jKatgert et al. 



2008. 



2001 



2009|). 



This has enabled to be tter understand simple shear phenomena, like localiza- 



tion or shear-banding (Wan g et al 



2008 



2006 



Janiaud et al. 



2006; 



Katgert et al. 



Langlois et all. 120081). alt h ough there are still open issues, such as the 



role of diso rder (IKatgert et al. 



2008, l2009h and that of bubble-bubble in- 



teractions (IDenkov et all |2009i ). Moreover, bubble monolayers are usually 



confined by solid plates, except bubble rafts, and the viscous friction between 



2 



bubbles and walls can drastically affect the flow profiles fjWang et all 120061 ) . 
This specificity of quasi-2D foams is a limitation for direct comparison with 
3D foam flows. 



Recently, tensorial co nstitutive mode 
coela s toplastic materia ls (ISaramito 



2008 



2007 



s for foams, or more generally vis- 



Benito et al. 



2008 



Cheddadi et al. 



Saramitd . 120091 ) . have been proposed in order t o ext end predictions 



of foam flows beyond pure shear (IMarmottant et all 120081 ) . They essen- 
tially rely on a local coupling between the local magnitude and orienta- 
tion of Tls, the deformation rate, and the elastic strain or stress. There- 
fore, to assess the validity of these tensorial models, local measurements of 
the elasticity, plasticity and flow profiles of foams flowing in complex ge- 
ometries is required. However, such flows have been less investigated than 
simple rheometric fl ows: some experiments h ave focuse d on flows past ob 



stacle s either in 3D (ICantat and Pitoisl . 120061 ) or in 2D (IDollet et al. 



2005 



20061). and some studie s have repor t ed on foam flows t hrough contractions 



flEarns 
2006). 



law and Wilsonl. 



1995, 



Asipauskas et al. 



200 



Bertho et al. 



1996 

Earnshaw and Wilsonl (jl995l . Il996l ) have shown that the rate of plas 



tic events is globally proporti onal to the local rate of strain, consistently 
with pure shear configurations (IGopal and Durianl . Il995l ). but in the absence 
of in-situ measurements of bubble deformation, the role of the elastic stress 
remains unclear. Hence, the existing studies have not provided a complete 
enough description of the flow to constitute stringent tests for the models. 

To help filling this gap, in this paper, a characterization of the rheological 
response of a quasi-2D foam flowing towards, through, and out a contraction 
is proposed. We introduce the experimental setup and the methods of image 
analysis in Sec. [2] We then describe in details a "reference" experiment, 
quantifying the elastic and plastic behaviors of the foam, and its flow profile 
(Sec. [3]). Several control parameters are then investigated, such as the flow 
rate, bubble size, contraction geometry and interfacial properties of the used 
surfactants (Sec. H|). We finally discuss and explain the complex interplays 
between elasticity, plasticity and flow, and between surface and bulk rheolo- 
gies (Sec. |5]), which make this configuration suited to shed insight on foam 
rheology. 
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2 Materials and methods 



2.1 Experimental setup 



We have used the foam channel fully described in lCantat et al.l (120061 ). It is 
a Hele-Shaw cell, made of two horizontal glass plates of length 170 cm and 
width 32 cm, separated by a gap h = 2 mm thin enough that the foam is 
confined as a bubble monolayer (Fig. [JJi). To make the contraction, home- 
made plastic plates were designed. They are long, narrow plates with a 
central wider part, which constitutes the contraction itself (Fig. [Tb). Three 
pairs were made, with different contraction lengths l c (2, 5 and 15 cm; the 
corresponding lengths L are 35, 33.5 and 28.5 cm). These plates are inserted 
in the channel through its open end, such that the distance between the 
contraction and the channel exits is 30 cm. This enables to change easily 
the contraction width w. Its maximal value is 4.4 cm when the plates are in 
contact with the channel side walls, and we studied also widths of 3.2, 2.1 
and 1.0 cm, for which a small gap remains between the channel side walls 
and the plates (Fig. [Tb). The corresponding aspect ratios W/w, with W the 
channel width upstream and downstream the contraction, are 4.6, 6.0, 8.6 
and 19. The dimensions and the alignment of the two plates of a given pair 
are within 1.5 mm, as can be seen in Fig. [Tb. 

The channel is connected upstream to a vertical chamber (Fig. []Ji) in 
which a given amount of soap solution is introduced thanks to a peristaltic 
pump. Nitrogen is continuously blown through injectors at the bottom of this 
chamber, producing rather monodisperse bubbles (Fig. [Tb) in the flow rate 
operating range of the paper (less than 150 ml/min per injector). The flow 
rate in each injector is independently controlled with an electronic flow-rate 
controller (Brooks). The resulting foam accumulates on top of the chamber, 
over a vertical distance where it drains, then is pushed through the channel. 
The transit time through the whole channel is less than 10 minutes in all 
experiments; we did not observe significant change of bubble size during 
this time, hence coarsening is negligible. The level of the foam/solution 
interface in the chamber was kept constant (within 2 mm) in all experiments 
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Figure 1: (a) Side view of the setup (without the contraction plates), (b) 
Top view (not to scale) of the foam channel and the contraction plates. Only 
the invariant dimensions have been indicated. The streamwise axis x, and 
the spanwise axis y, originating from the contraction center, are also drawn 
in dashed lines, (c) Snapshot of a foam flowing through a contraction of 
length 2.0 cm and width 3.2 cm. (d) Skeletonized version of snapshot (c), 
used for image analysis. The dashed rectangle is one of the boxes over which 
velocity, elastic stresses and plastic events are counted (see Sec. 12.21) . The 
flow direction is depicted by the bold arrows in each subfigure. 
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to minimize the variations of liquid fraction. The latter is very low (Fig. [Tb) 
and difficult to measure with precision. We can only give a rough estimate 
based on the decrease of the level of the foam/solution interface, which gives 
the volume of solution evacuated from the vertical chamber within the flowing 
foam. Comparing this volume with the gas flow rate, we get a liquid fraction 
between 0.2% and 0.4% for the experiments presented in this paper. 

In order to change the surface rheology, two different solutions were used. 
Most of the experiments were done with a solution of SDS (Sigma- Aldrich; 
purity > 99.0%) dissolved in ultra-pure water (Millipore) at a concentra- 
tion of 10 g/1, above the critical micellar concentration (cmc) of 2.3 g/1. 
Its surface static and dynamic properties were measured with a tensiometer 
(Teclis). We measured by the rising bubble method a the surface tension 
7 = 36.8 ± 0.3 mN/m. Using the oscillating bubble method, we measured a 
total dilatational surface modulus, defined as E$ = |d7/dlnA| with A the 
bubble area, below the noise level, hence lower than 1 mN/m: as expected, 



pure SDS interfac es have a neg 



as fully "mobile" (IDenkov et al. 



igible viscoelasticity and can be considered 



20051 ). All SDS solutions were used within 



a day from fabrication. Conversely, to study the case highly viscoelastic 

SLES, CAPB and myr istic acid 



Golemanov et al. 



(J2008|): we pre- 



or "rigid" interfaces, we used a mixture o: 
(MAc) , following the protocol described in 
pare a concentrated solution of 6.6% wt of SLES and 3.4% of CAPB in 
ultra-pure water, we dissolve 0.4% wt of MAc by continuously stirring and 
heating at 60° C for about one hour, and we dilute 20 times in ultra-pure 
water. The solution has a surface tension of 22.0 ± 0.5 mN/m, and a surface 
modulus of 216 mN/m for a frequency of 0.2 Hz and a relative area varia- 



tion 5S/S = 1.0%. Althou g h it 
"rigid" limit J 
measured by 



Denkov et al 



ully ensures that the interfaces are in the 
2005). this value is somewhat lower than that 



Golemanov et al 



(12008h . which is E s = 305 mN/m. This likely 
comes from the fact that these authors have used a piezoelectric control of 
the bubble oscillations, which was shown to be of better accuracy to impose 
sinusoidal area variations than o ur syringe- driven co ntrol, especially for in 



Golemanov et al. 



terfaces of high surface modulus fjRussev et all 120081 ). As 
(120081 ). we measured a significant decrease of the surface modulus with in- 
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creasing relative area variation (data not shown). Both solutions have a bulk 
viscosity equal to that of water, /i = 10~ 3 Pa s. 

The contraction region is lit by a circular neon tube of diameter 40 cm, 
placed just below the channel on a black board. It gives an isotropic and 
nearly homogeneous illumination over a diameter of about 20 cm. Movies of 
the foam flow are recorded with the high-speed camera APX-RS (Photron) 
at a frame rate of 60 or 125 frames per second, with a short shutter speed of 
1 ms, so that even the fastest bubbles (more than 10 cm/s) remain sharp. For 
all experiments except the ones with varying width (Sec. 14.21) and surfactants 
(Sec. I4.3p . we have recorded two movies, one upstream and one downstream 
the contraction, to record the flow far enough the contraction. The movies 
are constituted of 1000 images constituted by 1024 x 1024 pixels. 



2.2 Image analysis 



To extract the relevant rheological informa tion from the movies, we follow a 
procedure very similar to that presented in lDollet and Granerl (120071 ). First, 
each image of the movie is thresholded and skeletonized with a custom Im- 
ageJ macro. Since the foam is very dry and therefore the liquid films are 



thin, t he shape of the bu 



foams f Dollet and Graner 



Dies is well preserved (Fig. [TH), contrary to wet 



20071 ). As an important consequence, the elastic 



stress, based on the network of bubble edges, can be precisely estimated. 
The bubbles touching the contraction plates (sketched in white in Fig. [Tb) 
are less well preserved, because of the uncertainty on the location of these 
plates; notably, their area is often underestimated. Therefore, we will dis- 
card most of the information close to the contraction plates. Second, the 
sk eletonized movie is analyz ed by a custom Delphi program fully described 



in 



Dollet and Granerl (120071 ); based on individual bubble, edge and vertex 



tracking, it enables to compute the velocity, elastic stress and plastic events 
over a rectangular mesh of 30 x 30 boxes covering each image. The area of 
each box is comparable with the bubble area (Fig. [T]l), but the average over 
the 1000 images of each movie ensures that the fields are smooth at such a 
length scale. For the sake of clarity, the maps of the various fields presented 
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in this paper are recomputed and displayed at a coarser scale. 

The velocity is computed by averaging every bubble displacement over 
consecutive images. We will plot it both as a vector field v, and as a stream- 
line plot. We also compute the deformation rate tensor D = (Vi> + *Vt>)/2, 
where * designs matrix transpose. 



To compute the elastic stress flBatchelorl. Il970l). w e use the specific defi- 
nition valid for 2D foams ( Janiaud and Graner . 12005 ): 



a = \p 



with A the effective line tension (i.e. the pulling force exerted by each bubble 
edge), p the areal bubble edge density, i the notation for the vector joining 
two neighboring vertices (the edge curvature is neglected in the above ex- 
pression of the elastic stress), and (g) the tensorial product: in components 
a, j3 — x, y, Eq. (TTj) can be rewritten as a a p = \p{l a tpl t) . Each bubble edge 
is constituted of a thin liquid film separated by two parallel vertical inter- 
faces, between two Plateau borders in contact with the top and bottom walls. 
Since the foam is dry, the Plateau borders are of much smaller size than the 
gap, so we take the approximation: A ~ 27/i. Concerning the edge density, 
since in average a bubble has six edges, we take p = 3/ A with A the aver- 
age bubble area. The elastic stress, a purely mechanical notion, is strongly 
correlated to the p urely geometri c al not ion of bubble (elastic) deformation, 



cases (lAsipauskas et al. 



200 



3; 



as wa s proposed in | Aubouv et al.l fl2003[) and shown in various experimental 



Janiaud and Graner 



2005 



Marmottant et all 



20081 ): hence, we may hereafter interpret some features of the elastic stress in 
terms of bubble deformation. By definition, the elastic stress is a symmetric 
tensor with positive eigenvalues; hence, it is represented as an ellipse which 
major (minor) axis is proportional to the highest (lowest) eigenvalues and 

oriented along the corresponding eigenve ctors. 

The Tls are tracked as described in iDollet and Granerl (120071 ). For the 



four bubbles concerned by a Tl, we denote r<± the vector linking the centers 
of the two bubbles that lose contact, and f a the vector linking the centers of 



S 



the two bubbles that come into contact, and we ascribe this information to 
the box where the event takes place. In our routine, the research of appearing 
and disappearing contacts are run independently; therefore, to ensure that 
it is a relevant way to characterize Tls, we have checked that the number of 
disappearing (Nd) and appearing (N a ) edges is equal (within 10%) in each 
box. We thus compute the scalar field of the frequency of Tls per unit time 
and area: 

, N a + N d 

JTI 



94, i ■ 

i,x:1 box''movie 

where A^ ox is the area of a box and i mov ie the duration of a movie. 

As a limitation of our approach, we cannot measure the pressure field. 
Laplace law, which states that the pressure differences between two neighbor- 
ing bubbles is proportional to the curvature of their common edge, cannot be 
used here: it requires to know the full 3D geometry of the bubbles, whereas 
the out-of-plane edge curvature is lost during our purely 2D image analysis. 
Furthermore, Laplace law is strictly valid only at equilibrium, and its appli- 
cability to flow situations is questionable. Nevertheless, we verified in the 
reference experiment that the flow rate across various cross-sections remains 
constant (within 4%), which is a signature of an incompressible flow. Since 
the trace of the deformation rate tensor, V • v, is zero for such a flow, the 
deformation rate is a symmetric and (almost) traceless tensor, hence it has 
a positive and a negative eigenvalue: we represent two orthogonal lines, a 
thick (thin) one in the direction of the eigenvector associated to the positive 
(negative) eigenvalue, i.e in the dir ection of local elongation (compression) 



rate, as in 



Dollet and Grane 



3 koO^ \. 



3 Study of a reference experiment 

We now study in detail a reference experiment. It is a flow of a SDS foam, at 
flow rate 150 ml/min, in a contraction of length 2 cm and width 3.2 cm. The 
average bubble area is 39 mm 2 , and the polydispersity index, defined as the 
standard deviation of the list of individual bubble areas, is 22%, but 90% of 
the bubbles are within 3% from the average area, the rest being constituted 
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mainly of much smaller bubbles. Hence, the foam is rather monodisperse 

(Fig. db). 

We describe the main features of each field (velocity, Sec. I3.1| elastic 
stress, Sec. 13.21 and Tls, Sec. 13. 3p on the maps, and we report the evolution 
of the different components of each field along different axes: the central axis 
y = 0; an off-centered, spanwise axis located half-way between the central 
axis and the side walls, i.e. at \y\ = 4.8 cm, with an average over the axes 
y = —4.8; cm and y = 4.8 cm; a spanwise axis across the flow upstream the 
contraction, at x = —9.6 cm; and the symmetric spanwise axis downstream 
the contraction, at x — 9.6 cm. For the two latter axes, since y = is an 
axis of symmetry, the fields will be plotted for \y\ > 0, with an average over 
— y and y. 



3.1 Velocity 

The velocity field and the streamlines are plotted in Fig. [2] and the velocity 
components along the different axes are plotted in Fig. [3j The foam tends 
towards a plug flow far from the contraction, and is faster close to the con- 
traction, as can be seen from the velocity variation along the central axis 
(Fig. [3^). We do not evidence any dead zones, nor vortices, in the corners 
at the entrance side (denoted by letters A and B in Fig. [lb); the bubbles 
keep flowing towards the contraction even in the corners. This is a qualita- 
tive difference with pure viscoplastic fluids, w here dead zones or vortices are 



observed in the corners at the entrance side rtAbdali et all Il992l ; iJav et al. 



2002j). Moreover, several features show that the flow displays a strong fore- 



aft asymmetry with respect to the contraction center. First, the flow along 
the channel walls y = 1 cm is faster than along the channel walls y = —1 cm. 
Second, whereas the flow is purely convergent towards the contraction at the 
entrance zone, i.e. v y has always the opposite sign as y (Fig. [3b), it is not 
purely divergent at the exit: Fig. [2] shows that any streamline passes by a 
maximum distance from the central axis, before converging again towards 
the center. This reversal of the spanwise velocity component is clearly seen 
in both on an off-centered streamwise axis (Fig. [3b) and across a spanwise 
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Figure 2: Velocity field and streamlines for the reference experiment. The 
points where the velocity could not be reliably evaluated are left in blank. 
The streamlines starting close to the contraction plates are interrupted in 
zones where the velocity could not be reliably evaluated. 

axis in the exit region (Fig. [3H). Third, whereas the streamwise component 
of velocity along the central axis continuously increases towards the contrac- 
tion in the entrance region, it passes through an undershoot (at x — 10 cm) 
in the exit region (Fig. [3^); at this point, v x is 20% lower than the plug 
flow velocity. Furthermore, whereas the streamwise velocity component de- 
creases from the central axis to the sides in the entrance region (Fig. [3fc), 
it increases in the exit region and almost doubles from the central axis to 
the sides (Fig. [3H). We finally notice that along the central axis, v y slightly 
deviates from the zero value expected from symmetry. This is likely due the 
small geometric imperfections of the contraction (Sec. 12.11 and Fig. [It). 

3.2 Elastic stress 

The map of elastic stress is presented in Fig. H^,. Not surprisingly, the maxi- 
mal stress is directed towards the contraction at the entrance region, because 
the bubbles tend to be deformed by the converging flow. The situation is 
more complex at the exit region: the bubbles relax and revert their elas- 
tic stress very abruptly just at the contraction exit, and they are strongly 
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Figure 3: Plots of the velocity components v x (plain curve) and v y (dashed 
curve) along (a) the central axis y = 0, (b) the off-centered axis y = 4.8 cm, 
(c) upstream the contraction, along the axis y = —9.6 cm, and (d) down- 
stream the contraction, along the axis y = 9.6 cm. 



stressed along the orthoradial direction (from the contraction center) up to 
a few centimeters downstream the contraction, before experiencing a grad- 
ual elastic relaxation towards equilibrium, which is not finished as they are 
advected away from the observation window. As a specific feature of the 
setup, Fig. Hb shows that the foam is prestressed at the upstream end of 
the observation window; indeed, the foam has passed a first contraction from 
the main foam channel (of width 32 cm) to the channel of width W (Fig. 
[[]) 15 cm only before the upstream end of the observation window, hence it 
probably did not relax completely before arriving in the field of view. 

3.3 Tls 

We now turn to the plastic rearrangements Tls. The spatial distribution 
of their frequency is plotted in Fig. [5^- It shows that plasticity occurs 
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Figure 4: (a) Elastic stress field for the reference experiment. The major (mi- 
nor) axis represents the local direction and magnitude of maximal (minimal) 
stress, see Sec. I2.2I The points where the elastic stress could not be reliably 
evaluated are left in blank. Plots of the elastic stress components a xx (plain 
curve), o xy (dotted curve) and a yy (dashed curve) along (b) the central axis 
y = 0, (c) the off-centered axis y = 4.8 cm, (d) upstream the contraction, 
along the axis y = —9.6 cm, and (e) downstream the contraction, along the 
axis y = 9.6 cm. 



rather close to the contraction, within a radius of 7 cm. This is confirmed by 
Figs. Eb, c, d and e: there are about 5 times less Tls along the off-centered 
streamwise axis than along the central axis, and more than 10 times less 
across the spanwise axes located 9.6 cm from the contraction center. As an 
immediate consequence, the foam behaves as a viscoelastic medium outside 
this zone, which confirms that the features reported farther downstream (re- 
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focussing of the streamlines towards the central axis, velocity undershoot) are 
likely due to the elastic nature of the foam. Fig. [5^ and b show that there 
are more Tls at the entrance than at the exit of the contraction; however, 
this may be a specific result of the setup due to the prestress of the foam 
arriving in the contraction, which facilitates the occurrence of Tls. The 
distribution of Tls displays more generic fore-aft asymmetries: there are 
significant secondary maxima close to the corners in the exit region, and not 
in those of the entrance region; and, remarkably, there is a small zone of 
negligible plasticity just at the exit of the contraction (Fig. [5^ and b), which 
we will describe further in Sec. [3J 

4 Influence of control parameters 

We now study the influence of several control parameters on the flow of foam 
through a contraction. Starting from the reference experiment studied in 
Sec. El we show that a change in flow rate and in bubble area does not 
modify the velocity and elastic stress, up to a rescaling by the flow rate and 
bubble area, respectively (Sec. 14. lj) . We then study the influence of the 
length and width of the contraction (Sec. \4.2\i . We finally show how the 
physico-chemistry, i.e. the used surfactants, can affect the foam response in 
velocity and elastic stress (Sec. I4.3p . A summary of the various experiments, 
and their main characteristics: applied flow rate, maximal velocity, soap 
solution, mean bubble area, polydispersity, contraction width w and length 
£ c , is displayed in Tab. [H 

4.1 Flow rate and bubble area: flow rescaling 

To study the influence of the applied flow rate starting from the reference 
experiment of Sec. [31 for which the gas was injected through a single injector 
at a flow rate of 150 ml/min, we used a second injector through which gas is 
injected at the same flow rate, independently controlled by a second flow-rate 
controller. To study the influence of the bubble area at the reference flow 
rate of 150 ml/min, since the flow rate per injector is the key parameter to 
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Figure 5: (a) Spatial distribution of the frequency of Tls, fci, expressed 
in cm _2 -s _1 , for the reference experiments. Plots of the frequency of Tls 
along (b) the central axis y — 0, (c) the off-centered axis y = 4.8 cm, (d) 
upstream the contraction, along the axis y = —9.6 cm, and (e) downstream 
the contraction, along the axis y = 9.6 cm. 



tune the bubble size, we blow gas through the two injectors at 75 ml/min 
each. We also performed an experiment with a single injector blowing at 
75 ml/min. Note that we were limited to a narrow range of parameters by 
the number of injectors bubbling identically. 

To compare these four experiments, we focus on the variations on the 
central axis of the relevant components of the fields, properly rescaled. For 
the velocity, we take v x /vo, where v$ is the entrance velocity, i.e. the velocity 
averaged over a cross-stream section at the upstream end of the field of view. 
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4.2 


150 


3.90 


SDS 


33 


15 


5.0 


3.2 


4.2 


3 x 150 


9.84 


SDS 


34 


36 


15.0 


3.2 


4.3 


20 


0.26 


* 


19 


36 


2.0 


3.2 


4.3 


40 


0.58 


* 


19 


36 


2.0 


3.2 


4.3 


70 


1.02 


* 


19 


36 


2.0 


3.2 



Table 1: Summary of the main characteristics of the experiments presented 
in this paper. The star designs the SLES/CAPB/MAc solution described 
in Sec. 12. 1[ The fluid fraction is comprised between 0.2% and 0.4% in all 
experiments (Sec. 12. ip . 



For the 75, 150, 2 x 75 and 2 x 150 experiments, we found respectively 
entrance velocities of 0.29, 0.58, 0.62 and 1.15 cm/s. For the elastic stress, 
the relevant components are the normal components a xx and a yy . We will 
consider the normal stress difference a xx — a yy . Since the elastic stress scales 
as a bubble characteristic length, we rescale the normal stress difference by 
the trace a xx + a yy , to compare different bubble areas. 

We thus plot along the central axis v x /v , [a xx — a yy ) / [a xx + a yy ) and the 
Tl frequency for the four experiments, in Fig. [6j Fig. [6^ shows that all data 
for the rescaled velocity collapse on a single master curve; notably, the veloc- 
ity undershoot at x — 10 cm is a robust observation. Fig. [6]d shows that the 
rescaling is good also for the dimensionless normal stress difference, within 
slightly larger, but apparently random, deviations from the main trend. Since 
the elastic stress depends on individual bubble sizes, this may be due to the 
small variations of polydispersity between the different runs. However, we 
consistently observe the same variations of the dimensionless normal stress 
difference along the central axis: a small increase, then a plateau close to 
the contraction entrance, followed by a quick reversal just at the contraction 
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x (cm) x (cm) x (cm) 



Figure 6: Comparison between different flow rates and bubble areas: in terms 
of number of injectors times the flow rate per injector (which specifies see 
bubble area, see text for details) 150 ml/min (reference experiment, plain 
curve), 2 x 150 ml/min (dashed curve), 2 x 75 ml/min (dotted curve), and 
75 ml/min (dash-dotted curve). Plots along the central axis x = 0, of (a) the 
streamwise velocity component rescaled by the entrance velocity (see text for 
details), v x /vo, (b) the normal stress difference rescaled by the total elastic 
stress, (cr xx — CF yy )/{a xx + cr yy ), and (c) the frequency of Tls. 

exit followed by an elastic relaxation. Finally, the Tl frequency shows con- 
sistently two peaks, one higher at the entrance and one smaller at the exit, 
separated by a narrow zone with negligible plasticity just at the exit of the 
contraction (Fig. [6b). The Tl frequency for the 2 x 75 and 2 x 150 experi- 
ments is about the same, and is then lower for the 150 experiment, and even 
lower for the 75 experiment. This suggests that the Tl frequency increases 
as the velocity increases and the bubble size decreases, and dimensionally, we 
could expect that fri oc v/A 3 ^ 2 . But if this relation holded, the ratio of the 
Tl frequency of the 2 x 150 and 2 x 75 experiments would equal 2, whereas 
Fig. [6b shows that they are almost equal. Hence, there is no obvious scaling 
for the Tl frequency. 

4.2 Geometric parameters 

We now describe the influence of the geometric parameters of the contraction 
itself, i.e. its width and length. First of all, starting from the reference 
experiment of Sec. [31 we have changed only the width, and studied four 
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values of it: 1.0, 2.1, 3.2 (reference case) and 4.4 cm. As in the reference 
experiment, we have used one injector at 150 ml/min. 

As in Sec. 14. 1\ a proper rescaling is useful to compare data on the veloc- 
ity and elastic stress. A reference velocity based on the common flow rate 
of 150 ml/min is irrelevant, since as the contraction width changes, there 
is a varying gap between the contraction plates and the channel side walls 
(Fig. [T]) hence a varying "leakage" flow rate through this gap. Moreover, for 
these experiments, we have recorded only one movie each, centered on the 
contraction, hence we resolve the flow on a narrower range up- and down- 
stream the contraction (|x| < 9 cm) and we cannot estimate an entrance 
velocity in a plug flow regime as in Sec. 14.11 Therefore, we take as reference 
velocity the maximal velocity (see Tab. [T] for values), which is reached at 
the contraction center for these four experiments. Since the average area 
varies between the different runs (Tab. [T]), we also rescale the stress as in 
Sec. 14.11 From the results of last Section, we expect differences between the 
four experiments to be purely due to the width variation. Hence, we plot 
v x /v(x = 0), (a xx — cr yy )/(a xx + a yy ) and the Tl frequency along the central 
axis for the four experiments, in Fig. [7J 

The modification of the width does not alter the qualitative trends of the 
fields. As before, the velocity passes by a maximum in the contraction (and 
we miss here the undershoot, too far away downstream), the stress reverses in 
the contraction, and the Tl frequency shows two peaks at the entrance and 
at the exit, the former being slightly higher than the latter. Basically, the 
smaller the width, the more abrupt the velocity and stress variations across 
the contraction, and the higher the maxima in the Tl frequency distribution, 
which is coherent since the deformation rate increases. 

Second, using different pairs of plates, we have varied the length of the 
contraction: from 2 cm (reference case), to 5 and 15 cm. The width is 
maintained at the reference value of 3.2 cm. As for the different widths, we 
choose to rescale the velocity by v(x = 0), equal to 3.97, 3.90 and 9.84 cm/s 
respectively for the lengths 2, 5 and 15 cm, and we plot v x /v(x = 0), (a xx — 
<J yy)/(°'xx + Cyy) and the Tl frequency along the central axis for the three 
experiments, in Fig. [SJ 
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Figure 7: Comparison between different contraction widths, at the reference 
flow rate and bubble area: 1.0 cm (dashed curve), 2.1 cm (dotted curve), 
3.2 cm (reference case, plain curve) and 4.4 cm (dash-dotted curve). Plots 
along the central axis x — 0, of (a) the streamwise velocity component 
rescaled by the maximal velocity, v x /v(x = 0), (b) the normal stress dif- 
ference rescaled by the total elastic stress, (a xx — a yy )/(a xx + cr yy ), and (c) 
the frequency of Tls. 

Fig. [H^ shows a new interesting feature: for long enough contractions, 
the velocity is not maximum anymore at the center of the contraction, as 
was the case up to now for the "short" contraction of 2 cm. Namely, for a 
contraction length of 5 cm, the velocity is maximum at x = 1.3 cm, close 
to the contraction exit. For a contraction length of 15 cm, two velocity 
overshoots appears, at the entrance (x = —6.0 cm) and exit (x = 6.4 cm) of 
the contraction; these two overshoots are respectively 7% and 9% higher than 
v (x = 0). Let us notice also that the velocity inside the contraction does not 
reach a plateau, even in the longest contraction, hence its length/width ratio 
(equal to 4.7) is not high enough for a plug flow to be fully established. 

Concerning the stress, Fig. [Sb shows that there is a slow relaxation in the 
two longest contractions, but not big enough to come back to equilibrium. 
For the contraction length of 15 cm, there is also a small overshoot of (a xx — 
a yy)/{ a xx + &yy) before its fast reversal, which remains a very robust feature 
of all the studied flows. The two maxima for the Tls, one stronger at the 
entrance and one slightly lower at the exit, remain for the contraction lengths 
of 5 and 15 cm; moreover, there is in both cases a bit of plasticity also within 
the contraction, associated to the small stress relaxation. 
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Figure 8: Comparison between different contraction lengths, at the reference 
flow rate, bubble area and width: 2 cm (reference case, plain curve), 5 cm 
(dashed curve) and 15 cm (dotted curve). Plots along the central axis x = 0, 
of (a) the streamwise velocity component rescaled by the velocity at the mid- 
dle of the contraction, v x /v(x = 0), (b) the normal stress difference rescaled 
by the total elastic stress, (a xx — cr yy )/(a xx + cr yy ), and (c) the frequency of 
Tls. 



4.3 Surfactants 

Up to now, all presented experiments have been performed with a SDS so- 
lution, with a negligible dilatational surface viscoelasticity. We now study 
a solution of the surfactant mixture SLES/CAPB/MAc (see Sec. 12. ip . with 
a high dilatational surface viscoelasticity, to quantify the interplay between 
surface and bulk rheologies. The geometry of the contraction is that of the 
reference experiment: length 2 cm, and width 3.2 cm. Since we expect the 
velocity to be the key parameter, we have tried to vary the applied flow rate 
at given bubble area. To do so, we have prepared a foam in the whole chan- 
nel at a given flow rate of 70 ml/min, and we have changed the flow rate 
just before recording; thanks to the big length of the whole foam channel, 
the recorded flow is still made of bubbles generated with the flow rate of 
70 ml/min. We have performed three experiments, with maximal velocity 
(always reached at the contraction center): 0.26, 0.58 and 1.02 cm/s. We 
could not reach higher velocities, because the foam then collapses according 
to a process which will be described in subsequent studies. Since we want to 
track small structural changes, we have zoomed into the contraction, to get a 
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better resolution of the bubble shape at the expense of a smaller field of view 
(about 7x7 cm 2 centered on the contraction). In order to compare these 
experiments with the reference experiment of Sec. El we plot the streamwise 
component of velocity rescaled by the maximal velocity, v x /v(x = 0), and 
the dimensionless stress difference (a xx — a yy )/(a xx + cr yy ), along the central 
axis. In the narrow window of observation considered here, the statistics of 
plastic events becomes too poor to be studied here. 

Fig. [9^ shows that the rescaled velocity does not vary much between the 
different experiments, and that the deviations from the reference experiment 
remains small. On the contrary, Fig. [Hb shows a strong and systematic 
behavior difference between the SDS solution and the SLES/CAPB/MAc 
mixture: with the latter, higher values of the dimensionless stress difference 
{pxx — &yy)/{.&xx + ffyy) are reached, and the variation amplitude of this 
parameter increases with increasing velocity. As stated in Sec. 12.21 this 
is correlated with a bigger bubble deformation, which is readily seen by 
comparison between Figs. [9fc and d. 



5 Discussion 

5.1 Interplay between elasticity, plasticity and flow 

In Sec. [3] and HI we have described separately the behavior of the veloc- 
ity, elastic stress and Tl fields. We now discuss the coupling between these 
quantities. First, the modification of various control parameters, notably 
flow rate, bubble area and contraction geometry (Sec. 14.11 and |4.2[) . has been 
shown not to change the qualitative features of the flow of foam through a 
contraction. The effect of the foam disorder on the flow pattern remains an 



open question, as was pointed out by lKatgert et al.l (120081 ). We have worked 



with rather monodisperse foams (Fig. [Tb), but our polydispersity indices 
are not negligible, and they vary significantly between various experiments 
(Tab. [T]). Hence, the qualitative trends of the flow of foam through a con- 
traction are not modified when a moderate amount of disorder is present. 
Consequently, the main features of the foam flow are direct consequences 
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Figure 9: Effect of the surfactant. Plot of (a) the streamwise veloc- 
ity component rescaled by the velocity at the middle of the contraction, 
v x /v(x = 0), and (b) the normal stress difference rescaled by the total elastic 
stress, [a xx — (J y y)/(<J X x + Oyy)-, f° r f° ur experiments: the reference experi- 
ment with a SDS solution (dotted curve), and three experiments with the 
SLES/CAPB/MAc mixture, with maximal velocity 0.26 (dashed curve), 0.58 
(dash-dotted curve) and 1.02 cm/s (plain curve). To illustrate the change of 
bubble deformation, snapshots of (c) the reference experiments and (d) the 
fastest experiment with the SLES/CAPB/MAc mixture, taken just at the 
entrance of the contraction (—3 cm < x < —1.5 cm), are displayed. 
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Figure 10: Deformation rate for the reference experiment. The thick (thin) 
lines represent the local maximal elongation (compression) rate, see Sec. 12. 21 



of the interplay of elasticity, plasticity and flow, that we can study on the 
reference experiment without m uch loss of generalit y. 

Ex periments in pure shear flGopal and Durianl . Il995l ) and contraction 
flows (jEarnshaw and Wilson! . Il995l ) have shown that the rate of plasticity 
is globally proportional to the rate of strain. Contrary to these studies, we 
have access to the local values of the deformation rate, which is displayed in 
Fig. [TUJ From these local values, we can check whether we recover a propor- 
tionality between the rate of plasticity and the deformation rate, by plotting 
for each box the frequency of Tls versus the norm of the deformation rate 
tensor, |5| = sj {D : D)/2 = ^J(D 2 XX + 2D 2 xy + D 2 yy )/2, in Fig. [TTJ It shows 
that there is a proportionality between the rates of plasticity and deformation 
when these are not too large, but that there is a huge scatter at large rates 
of plasticity and deformation. Notably, several points at high deformation 
rate show very little plasticity; actually, they correspond to the small zone 
just at the contraction exit, where a sharp minimum of plasticity was already 
pinpointed in Sec. 13.31 Hence, the proportionality between rates of plasticity 
and deformation is not valid everywhere. 

We now test a more adva nced coupling between ela sticity, plasticity and 



deformation rate proposed by 



Marmottant et al. 



(120081 ) . These authors used 
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Figure 11: Frequency of Tls versus deformation rate. Each box gives one 
data point. The straight line is the best linear fit through all data. 



a plastic tensor defined as (IGraner et all 120081 ) : 



P=\fri{{{r d ® r d ) - (f a ® r a ) ) ■ M^ 1 + * [( (r d <g> r d ) - (f a <g> f a ) ) • M' 1 } }, (2) 



where M = (r <g> f) is the texture tensor (jAubouv et all 120031 ) . defined on 
the network of the vectors r linking centers of neighboring bubbles. The 
plastic tensor encompasses not only the frequency of Tls, but also there 
direction; namely, the positive (negative) eigenvalue of this (almost) traceless 
tensor is the preferential direction of bubble separation (attachment). As 
the deformation rate, we represent the plastic tensor by a pair of orthogonal 
lines, the thick (thin) one in the directi on of the positive (negativ e) eigenvalue 
(Fig. fT2l) . As a constitutive relation, iMarmottant et al.l (120081 ) conjectured 



that the plastic tensor P, the deviator of the elastic stress tensor, 



a d = a - - (ti a)I, 



(3) 
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Figure 12: Plastic tensor P for the reference experiment. The thick (thin) 
lines represent the direction of the eigenvector associated to the positive 
(negative) eigenvalues, and the line lengths are proportional to the absolute 
value of the eigenvalues. The points where the plastic tensor could not be 
reliably evaluated are left in blank. 



where I is the identity tensor, and the deformation rate D are related through: 



if D : a d < 
2\crd\ 



(4) 



where h, a function varying between and 1, would be the Heaviside function 
if the foam showed no plasticity below a yield stress oy and yielded abruptly 
at that yield stress. Notic e that Eq. (j3J) slightly differs from the original ver- 
sion of 



Marmottant et al. 



(120081 ). where an elastic strain U was used instead 



of a; but since the de viators an. and Ud have been shown to be proportional 



flAsipauskas et al. 



2003 



Janiaud and Graner 



2005 



Marmottant et al 



2008 



it is equivalent to use Ud or <r d in Eq. (J4]). 

The quantity D : a d is proportional to cos 29, with < 9 < n/2 the angle 
between the major axes of elastic stress and deformation rate; it is positive 
if 9 < 7r/4, e.g. if the two axes are aligned, and negative if 9 > 7r/4, e.g. if 
the two axes are proportional. If the two axes are aligned, D : = 2\D\ |<7d|, 
hence Eq. (jjl) predicts that the rate of plasticity should be proportional to the 
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deformation rate. Since in most regions of the flow the major axes of elastic 
stress and deformation rate are indeed well aligned (Fig. [10] and Hk), it likely 
explains why there is in average a good correlation between the frequency 
of Tls and the deformation rate. On the contrary, the relation implies 
that no plasticity should occur if the elastic stress and the deformation rate 
tend to act along opposite directions; Fig. [10] and UK show that this is the 
case just at the contraction exit, and indeed this corresponds to the area 
with a marked minimum of plasticity. In this zone, the spanwise elongation 
rate serves not to trigger plastic events, but to reverse the elastic stress from 
streamwise to spanwise. 

There is no model yet for the function h entering Eq. (j3J), hence it is 
difficult to test the full relation. On the other hand, it predicts that P should 
be aligned with ad rather than D, which we test: we compute, for each box, 
the angle between the major axes of the plastic tensor and elastic stress, 
6(P,a), and, for comparison, the corresponding angle for the plastic tensor 
and deformation rate, 9(P,D), and we plot the histograms of distribution 
of these two angles in Fig. [TS] It shows that the plastic tensor is indeed in 
average well aligned with the elastic stress than with the deformation rate; 
we have (9(P,a)) = 0.280 rad = 16°, and (9{P,D)) = 0.383 rad = 22°. This 
result is a strong indication that the plasticity rate is not only related to the 
deformation rate, but also to the elastic stress. We do not proceed here with 
the discussion of the interplay of elasticity, plasticity and flow and with tests 
of the models. The next step would be to compare the full elastic, plastic and 
velocity fields in our experiments with simulations and ah initio predictions 
of the flow, but this is beyond the scope of this paper and will be pursued in 
future studies. 

At about 7 cm away from the contraction, plasticity vanishes, hence the 
foam behaves as a viscoelastic medium. Therefore, the velocity undershoot 
(reached in a zone where almost no Tls occur, see Figs. [3^ and [5b) and the 
re-focussing of the streamlines appear as elastic effects. Ind eed, the velocity 
undershoot is reminiscent of the ne gative wake phenomenon (IHassagerl . 119791 1 
previously reported also in foams (IDollet and Granerl . 120071 ) , and attributed 
to elastic effects; at the location of the velocity undershoot, the bubbles are 
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Figure 13: (a) Histogram of the angle between the major axes of the plastic 
tensor and elastic stress, 9(P,a). (b) Histogram of the angle between the 
major axes of the plastic tensor and deformation rate, 9(P,D). Angles are 
expressed in radians and are defined in the range [0, 7r/2]. 



very elongated spanwise, and they then elastically recoil streamwise (Figs. 

and b). Similarly, the re-focussing of streamlines may be attributed to an 
elastic recoil of the bubbles. 

On the other hand, the velocity overshoots observed at the entrance and 
exit of long enough contractions (Fig. [8^l) appear in yielded regions with 
a strong plasticity rate (Fig. [Bb), hence they do not originate purely from 
foam elasticity. This is confirmed by the recent first observation of a ve- 
locity overshoot at the entrance of a contraction for a viscoplastic fluid, by 
Rabideau et all There does not seem to be obvious qualitative explanations 



for these velocity overshoots. 



5.2 Interplay between surface and bulk rheologies 

We have shown in Sec. 14.31 that the elastic stress depends significantly on 
the velocity, in the case of the SLES/CAPB/MAc mixture. This interest- 
ing observation shows that when the friction in liquid films, either between 
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neighboring bubbles or between bubbles and walls, becomes strong enough, 
a dynamic parameter (here, the local value of the velocity or of the deforma- 
tion rate) can influence the structure. This dynamically-induced deviation 
from a quasistatic foam structure (such as that fully characterized in Sec. [21 
which was shown in Sec. 14. II indeed not to depend on the flow rate) has been 
predicted theoretically by lKraynik and Hansen! (119871 ) and observed in other 



contexts, e.g. by the deviatio n from the equilibrium la ws governing the angles 
at which a film meets a wall (IDrenckhan et al.L 120051 ) . Actually, this also oc- 
curs here: Fig. [9H shows that some liquid films with the SLES / CAPB /MAc 
mixture seem thick, seen from above. This apparent increased thickness is 
nothing but the larger horizontal section of the films that exhibit a strong 
curvature across the gap between the two plates; this is also related to devi- 
ations from the equilibrium laws, because a liquid film between two bubbles 
at equilibrium is straight and perpendicular to the walls. 

It would be interesting then to compare our e xperimenta l snap shots with 
shape predictions from the viscous froth model (IKern et al.L I2004T). which is 



already used to study t he morphology of single films and bubbles (IGrassia et al. 



2008; 



Cox et al 



20091 ) under the effect of friction acting on Plateau borders, 



where a liquid film meets a wall. However, we have here two sources of fric- 
tion: between bubbles and walls, or between neighboring bubbles. When the 
latter dominates, which is usually the case in 3D foams, it has been observed 
that the yield strain, i.e. the applied strain at the onset of plast icity, is an 
increasing function of the applied strain rate (IRouver et al.L 120031 ) . This was 
indeed interpreted as a consequence of the internal viscous stresses. Work is 
in progress to study the influence of one or the other on dynamical modifi- 
cation of the foam structure, in simpler geometries. 

Finally, let us stress that such a variation of the elastic stress is a priori 
possible also for the SDS solution. However, its surface viscoelasticity is two 
orders of magnitude lower than that of the SLES / CAPB /MAc mixture, hence 
deviations from the quasistatic regime are not expected under velocities of 
order 1 m/s, which is beyond the range of flow rates that we can produce 
with our setup. 
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6 Conclusions 

We have experimentally investigated the two-dimensional flow of foam through 
a contraction. Thanks to image analysis, we could extract a rather complete 
information on the foam behavior, in terms of elastic stresses, plastic events, 
and flow. The flow exhibits a strong fore-aft asymmetry (with respect to the 
contraction center) and some striking features, such as a velocity undershoot 
at the contraction exit. The qualitative features of the flow are rather insen- 
sitive to control parameters such as the bubble size and contraction width 
and length, as well as the flow rate, provided it is low enough (quasistatic 
regime). Conversely, as friction within the liquid films becomes significant, 
for instance using surfactants with a strong surface viscoelasticity, the elastic 
stresses show a strong dependence on the flow velocity. 

By comparing elastic and plastic local descriptors of the foam and its 
flow profile, we have been able to show that there is an essential coupling 
between elasticity, plasticity and flow. Notably, there is not everywhere a 
proportionality between the rates of plasticity and deformation, and the rel- 
ative orientation of elastic stresses and deformation rate plays a fundamental 
role, which encourages to further develop and test tensorial viscoelastoplastic 
models for foam rheology. As another perspective of this work, the inter- 
esting coupling between surface and bulk rheologies deserves to be further 
investigated and understood, first in simple geometries (plug flow). 
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